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We extract the K + K~~ couplings of the isoscalar scalar mesons <r//o(600) and /o(980) from -kit — » tttv/KK scatterings and 
found: \g a K+K- I/ISo-tt+tt- I — 0.8 and [p/gK+Jf- l/lfl/ ir+ir- I — I-?- These results, together with the tiny "direct" 77 width of 
the a and its large hadronic width, are a strong indication for the gluonium/glueball nature of the a-meson, as predicted by 
QCD spectral sum rules (QSSR) © some low-energy theorems (LET), while some other assignements (nn molecule, tetraquark 
state and ordinary qq meson) do not satisfy simultaneously these requirements from the data. These properties suggest that 
the a can be a scalar meson associated to the U(l)v conformal anomaly like is the jy'-meson for the U(1)a anomaly. 



1. Introduction 

Understanding the nature of scalar mesons in terms of 
quark and gluon constituents is a long standing puzzle 
in QCD [l|2j . The problem here is that some states are 
very broad (a and «(if confirmed) mesons) and others 
are close to an inelastic threshold (/o(980), an (980)), 
which makes their interpretation difficult. One might 
expect that the decay rate of these mesons into two pho- 
tons could provide an important information about their 
intrinsic composite structure. Indeed, a recent analysis 
of 77 — > 7T7r data [3] indicates that the cr/fo (600) (here- 
after called a) meson could be such a gluonic resonance. 
In this paper, we pursue the test of its nature by study- 
ing its couplings (and that of the /n(980)) to K + K~ 
and 7r + 7r~. 

2. Gluonium nature of the a from 7777/77 ~~ * 7r7r 
The existence of glueballs/gluonia is a characteristic pre- 
diction of QCD and some scenarios have been developed 
already back in 1975 @J. Today, there is agreement that 
such states exist in QCD and the lightest state has quan- 
tum numbers J PC = 0++. QSSR |5I6I7I8I9I2| determi- 
nations of its mass found 1 : 



M a 



1 GeV 



(1) 



confirmed recently by lattice simulation using dynami- 
cal fermions [10] and a strong coupling calculation [IT]. 
Some phenomenological implications of scalar glueball 
have been studied in the literature [51617181121131141151 
fTB"] . The analysis of the 77 — * 7777 data [3] using an 
improved model of [T7] leads to the (model-dependent) 
"partial" 77 widths: 



<t — >77 
presc 
a — >77 



(0.13 ±0.05) keV 
(2.7 ±0.4) keV , 



(2) 



where dir and resc refer respectively to the direct cou- 
pling of the a resonance and to the rescattering term ob- 
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1 QSSR also requires the existence of a higher gluonium mass of 
about (1.5-1.6) GeV from a consistency of the subtracted and 
unsubstractcd sum rules. 



tained using an unitarized Born amplitude. This leads 
to the (model-independent) total 77 width (direct + 
rescattering) : 



(3.9 ±0.6) keV , 



(3) 



which is in agreeement with the results from the ex- 
isting fits in [18I19I20I21I2"2"] . The previous results in 
Eqs. © and have been obtained at the complex 
pole obtained in [3] and given in Table [TJ with the cor- 
responding residue: 



0.06 - i 0.50 GeV 



(4) 



Table 1 

Mass and 1/2 width in McV of the a meson in the complex plane. 



Processes 



M a - iT a /2 



Refs. 



7777 — 


> tttt/KK 


422 - i 290 
441± 16 - i 272^5 
461 ± 15 -i (255 ±16) 


® 

m 


Jli>- 


-> CJ7777 


541 ± 39 - i (222 ±42) 


[25] 


D+ - 


■* 77 + 77 ~77 + 


478i^±17-i (162±^±21) 


26 



which is in agreement with the other results given there. 
These values of the a meson parameters in the complex 
plane, when appropriately translated into the real axis, 
become at the on-shell mass M° s w 0.92 GeV(see sec- 
tion [J): 



r° 



1.02 GeV 



-pos,dir 
a — »77 



(TO ±0.4) keV 



(5) 



which are in a remarkable agreement with the QSSR [27[ 
[28] and LET predictions obtained in the real axis 2 for a 
"bare/unmixed" gluonium/glueball ob state having the 
mass [5]: 



(0.95 - 1.10) GeV , 



(G) 



and the widths [61718] from the couplings in Table [U 

r„„^+„- ~ 0.5 GeV , 

~ (0.2 - 0.6) keV . (7) 



(7b — >7T+7T 

r 



2 These LET have been used earlier in |29I30I . 
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The large c-width into tttt indicates a strong viola- 
tion of the OZI rule in this channel and signals large 
non-perturbative effects in its treatment. This large 
hadronic width also disfavours its qq interpretation. 



Table 2 

QSSR ant LET predictions for the modulus of the hadronic couplings 
in GeV and 77 widths in keV of a scalar meson having a mass of 1 
GeV for gluonium I6I7I8J , qq meson I31I7I8I32I and four- quark states 



ED. 



Meson 


ffSV + 7I— 


gsK+K- 


r^77 


a B = gg 

S 2 = -^(uu- 


5 

1- dd) 2.5 


gaBK+Tr— 

1 

2yS , 27r + 7r~ 


0.2 - 0.6 

25 p 

9 L ao77 


S 3 = ss 


9S 3 K+K~ 


2.7±0.5 


±r 

g- 1 - 0077 

0.4 x 10~ 3 



In fact, using the hadronic couplings in Tabled QSSR 
predicts for a S2 = 1 / V2(uu + dd) , with mass of 1 GeV: 



\gs 2 ir+K- 1 

WttM^ 



4toJ 

02 



120 MeV 



(8) 



Using the o,q — > 77 width of about 2 keV from a QSSR 
analysis of a quark triangle loop vertex including the 
(ipip) condensate contribution [3 if 3 ! one predicts a S2 — * 
77 width of about 5 keV (an analogous result has been 
obtained in [55] Inland a S3 — > 77 width of about 0.2 keV 
(see also [35]). In the same way, QSSR also predicts, 
for a four-quark state, having the same mass of 1 GeV 
pM] , a 77 width of about 0.4 eV [3ljE These QSSR 
predictions and the value of the a — > 77 direct coupling 
from the data do not favour the and 4-quark scenar- 
ios. However, the conclusion is not sharp as some other 
approaches may still allow the possibility to have a four- 
quark state. 

For further tests of the nature of the a, we investi- 
gate the extraction of the a coupling to K + K~ . It is 
clear that, one expects a null value of this coupling in 
a tttt molecule and/or tetraquark assignements for the 
a. This is in contrast with its gluonium assignement, 
where one, instead, expects its large (almost) univer- 
sal coupling to pairs of pseudoscalar mesons |6|7|8j (see 
Table 



WaK+K- 



\g<r-K+-K- 



5 GeV 



(9) 



3. The g a K+K- coupling from tt + tt~ — * tttt/KK 
Unlike .9^+^-, this coupling cannot be directly mea- 
sured due to phase space suppression. In the following, 



3 A kaon tadpole loop mechanism leads to smaller value of the 
a -» 77 width of about 0.25 keV [12]. 

4 However, an approach based on the scalar anomaly leads to a 
smaller value of about 0.2 keV 1341 . 

5 An alternative approach based on kaon loop leads to a larger 
value of about (0.2 ~ 0.6) keV [37]. 

6 This typical non-pcrturbative prediction differs from a perturba- 
tive argument (which should not apply below 1 GeV) where the 
coupling behaves like the current quark mass. 



we extend the method in |17|3j used for elastic tttt scat- 
tering to determine the a parameters presented in the 
previous section. We shall also use S-matrix coupled- 
channel models with poles discussed in [39] . 



3.1. The analytic K-matrix model of [T7]l3] 

The strong processes are expressed by a K matrix model 
representing the tttt — > tttt/KK amplitudes by a set of 
resonance poles [IT]. In this case, the dispersion re- 
lations in the multi-channel case can be solved explic- 
itly, which is not possible otherwise. This model can 
be reproduced by a set of Feynman diagrams, including 
resonance (bare) couplings to tttt and K K and 4-point 
mr and KK interaction vertices. A subclass of bub- 
ble pion loop diagrams including resonance poles in the 
s-channel are resummed (unitarized Born). In [3], we 
study the elastic tttt scattering, where we introduce a 
shape function fo(s) = fp(s) which multiplies the otttt 
coupling and, for simplicity, we do not include the 4- 
point coupling term. Unlike approaches based on dis- 
persion relations, this approach can provide a separation 
of the "direct" resonance coupling with the "rescatter- 
ing" contributions which have been explicitly analyzed 
for 77 scattering in [3"pl . In the following, we discuss 
this approach (for a pedagocical reason) for the case of 
1 channel © 1 resonance. The real analytic function 
fp(s) : P = 7r, K is regular for s > and has a left cut 
for s < 0. For our low energy approach, a convenient 
approximation, which allows for a zero at s — sap and 
a pole at ap > simulating the left hand cut, is: 



sap 



fp{s) = 

S + (JDT) 

The unitary PP amplitude is then written as: 

T (0V S \ = Gpfp(s) = Gpfp(s) 
P s R -s-G P f P (s) V P {s) ' 



(10) 



(11) 



where the index corresponds to / 



0, 



e t5 p smS ( p'/pp(s) with p P (s) = (1 - Amp/s) 1 / 2 ; G P = 
gp B = <7^pp/(167r) are the bare coupling squared and 



Im V P = Im (-Gpfp) = -(9p P )Gp f P 



(12) 



with: (6pp)(s) — below and (9pp)(s) = pp(s) above 
threshold s — ' - 1 1 
where T>p(s 



T P °\s) 



9p 



4m P . The amplitude near the pole so 
and T> P (s) w 2? p (s )(s — s ) is: 

Gpfp(so) 



so 



9% 



(13) 



The real part of D P is obtained from a dispersion rela- 
tion with subtraction at s = and one obtains: 



fp{s) = ~ 

TT 



h P {s) -h P {0) 



(14) 



hp{s) = fp{s)L sl (s)-(a NP /(s+a D p))L sl (-app), <j np 
is the residue of fp(s) at —ctdp and: L sl (s) = 



7 A separation of the direct and rescattering term can also be stud- 
ied by measuring C assymetry in e+e — |38j . 
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[ (s - Amp) /m%]Li(s,mp) with L x from [T7]. This 
analysis can be generalized to the case of 2 channels 

2 resonance poles. A priori, the shape functions dif- 
fer for the 7T7T and K + K~ channels. However, evoking 
SU(3) symmetry, we can assume (to a first approxima- 
tion) that they are equal. 

3.2. Coupled channel model of |39j 

In the 2- and 3-coupled channel approach presented in 
[39] resonances correspond to the closest to physical re- 
gion poles of the S-matrix in complex energy plane. For 
each resonance such a pole is chosen among 2 % poles (i 
is number of coupled channels) what is in contrast with 
if -matrix models and those using Breit-Wigner formu- 
lae. In this way, the unitary S'-matrix approach in [39] 
delivers us a spectrum of scalar mesons below 1.6 GeV 
together with their couplings and branching ratios fitted 
to experimental data on the tttt and KK phase shifts 
and inelasticities. For example, for the Solution A in 
[39], the ratio of couplings to the KK and tttt channels 
for a and / (980) states were 0.25 and 2.24 respectively. 
These fits can be improved by using the method in [23] 
based on dispersive analysis of experimental data. The 
authors have shown that theoretical constraints given 
by Forward Dispersion Relations (FDR), sum rules and 
by once and twice subtracted dispersion relation (GKPY 
and the Roy's equations respectively) allow to determine 
the 7T7T scattering amplitudes consistent with analytic- 
ity, unitarity and crossing symmetry. 
Forward Dispersion Relations (calculated at t = 0) are 
used in [24] for three isospin combinations of tttt ampli- 
tudes: for 7r°7r + , 7r°7r° and for t-channel one with isospin 
It = 1. For example, the FDR for the former two com- 
binations which need two subtractions read: 
ReFi(s,0) - Fi(4ml,0) = 

OO 

s(s-iml) f (2s' - Aml)ImFi(s',0)ds' ( 15 ) 
w 7 s'(s' - s)(s' - 4m%)(s' + s - Aral) ' 

where Fi stands for 7r°7r + or 7r°7r° amplitudes, while 
FDR for the latter one do not need subtractions. In the 
ideal situation, the difference on the left hand side of 
Eq. (fTS"]) equals to zero but for realistic amplitudes this 
difference is minimized in the fitting procedure described 
in [24]. 

In the case of nonforward dispersion relation (the Roy's 
and GKPY equations) derived with imposed crossing 
symmetry condition, a difference between "output" and 
"input" partial wave amplitudes is minimized together 
with that for FDR. The "output" ones are calculated 
for three partial waves JI: SO, PI and S2 up to almost 

1 GeV and are given by 

Rc f I e (output \ S ) = ST(a° 0l a° 2 ) + 

Smax 

Jds'Kii:{s,s')imfi;{s')+di{s,s max ). ( 16 ) 

V V Ami 

The a[J and a° are the SO and S2 scattering lengths, 
KfJ (s, s') are kernels and d\ (s, s max ) are the so-called 



driving terms. The subtracting terms (ST) are linear 
combinations of scattering lengths in the GKPY equa- 
tions and do not depend on s. In the Roy's equations, 
these are also combinations of Oq and a® but depend lin- 
early on s — Am\. The integrals with kernels K\y (s, s') 
are calculated for partial waves with V < 2 for the Roy's 
equations and with I' < 4 for GKPY ones. The maximal 
value of the effective two pion mass squared s max , up 
to which experimental phase shifts and inelasticities for 
waves with I 1 and /' were parameterized, was chosen to 
be 1.42 GeV. The driving terms d\ (s, s max ) describe the 
influence of higher partial waves in the whole s range and 
of the lower ones above s = s ma x where Regge model 
was applied. The "input" amplitudes are those of which 
imaginary parts were used in once and twice subtracted 
dispersion relations in kernel and driving terms. 
As was shown in [40] the GKPY equations give output 
amplitudes with much smaller errors than those from the 
Roy's ones. It means that the GKPY equations impose 
stronger constraints on studied amplitudes. In addition 
to the Roy's and GKPY equations, two sum rules which 
relate high energy (Regge) parameters to low energy P 
and D waves were also considered. 
Following this method, we have performed similar fit for 
the SO amplitude from 2- and 3-coupled channel model 
of [39] , All other 7T7t partial waves were the same as in 
[24] and fixed. As a result, the ratio of couplings to the 
KK and tttt channels becomes 0.75 for a and 1.98 for 
/o(980), which we report in Tabled The change of this 
ratio for a in comparison with that for solution A of [39] 
confirms the general conclusion from [24] and [41j that 
available tttt experimental data sets, (also this from [42] 
used in fits of [39]) do not fulfill enough well theoreti- 
cal constraints from dispersion relations. It shows how 
important are theoretical demands given by dispersion 
relations and sum rules. Such constraints should be used 
in the fits of the experimental data of tttt amplitudes. 

3.3. Data input and results 

We have used the data of inelasticities and phase shifts 
from [43144145139] . More details of the analysis will be 
published elsewhere. In Table [3] we give the different 
couplings coming from the analyses based on the models 
in [17|3j and [39] . We shall use the normalization in Eq. 
p]). implying that |gp| 2 defined in section l3~Tl is equal to 
|5(tPp| 2 /(167t) in the narrow width approximation. We 
shall also use the relations: 

2 1 

\9Sir+ir-\ 2 = glSSTTTrl 2 , \gSK+K~? = ^SKK? -(17) 

We compare in the same Table [3] the results with the 
ones obtained in [22] from the use of dispersion relations 
in tttt and KK coupled-channel analysis in [IS] ■ We also 
compare the results with the ones from (j> — ► a/ / (980) 7 
[4"T] and J/ip — > 4> irir/KK [35] decays. One can notice 
that the different predictions for the a couplings are 
quite stable. The a couplings results using the model of 
[3] have been obtained in the case with 1 resonance © 
2 channels. In this case, one has obtained two output 
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poles but the 2nd one stays in an unphysical region. The 
results for the absolute value of the /o(980) couplings 
depend on the different models but the ratio of couplings 
is also quite stable. These results indicate that the a and 
/o(980) have important couplings to KK, where the one 
for the a is more remarkable. 



Table 3 

Modulus of the -k^tt and K complex couplings in GeV of the 
tr and of /q(980) from S- and K-matrix models for tttt — * tttt/KK 
scatterings compared with the ones from <f> and J/ift decays. rs^K = 
9 SK + K - /9stt+ tt— ' S = a, /q. means that this result has been 
obtained in the case 1 resonance © 2 channels. 



Processes 


9<T7r+TT 




r fa*K 


Models 


This work 












TTTT — * TTTT /KK 


1.55 


« 1 


1.57 


1.8 


[17|3| 




2.5 


0.75 




1.98 


m 


Others 












TTTT — > TTTT / KK 


2.5 


0.62 


1.55 


1.2 


m 


0->a// o (98O) 7 




0.67 






nu 


J/tp -> tttt/KK 






2.35 


1.8 


08] 


Average 


2.2 


0.8 


1.8 


1.7 





4. Comparison with the LET results 

4.1. LET for the a couplings to Goldstone bosons 

These couplings can be obtained from the vertex func- 
tion: 



V[q 2 = {q 1 -q 2 f] = {n 1 \e^ 2 ) 

where: 
1 



(18) 



E 

=u,d,s 



B^ = -P{a s )G^G%+ ^ (l+ 7 (a s ))m^ , (19) 



is the conformal anomaly (trace of the energy momen- 
tum tensor) with j3 and 7m are the QCD /3-function and 
mass-anomalous dimension while G^" is the gluon field 
strength and ipi is the quark field. V(q 2 ) obeys a once 
subtracted dispersion relation [BJ: 



V(q 2 ) = V(0) + q 2 



dt 



1 



4 m 'l 



t (t — q 2 — ie) tt 



-ImV(t) , (20) 



with the condition: V(0) — 0(m 2 ) 
limit. Using also the fact that V'(0) 
derive the two sum rules: 



b=<7B 



> in the chiral 
1, one can then 



where fs is the decay constant analogue to f v . The 
1st sum rule requires the existence of at least two res- 
onances coupled strongly to tttt. Considering the ctb 
and <j' b (its first radial excitation) but neglecting the 
small G-coupling to tttt as indicated by GAMS [49jj and 



8 The G(1.6) is expected to couple mainly to the U(1)a channels 
r/r/ and through mixing to 7/r/ , r\r\ |6I7| . 



some other data [50] , one predicts in the chiral limit [6|7j 
the almost universal couplings in Eq. © (see Tabled, 
which would correspond to a large width q 



0.5 GeV 



(22) 



This large width into tttt is a typical OZI-violation due 
to non-perturbative effects expected to be important in 
the region below 1 GeV, where perturbative arguments 
(like a vanishing of the hadronic glueball coupling in the 
chiral limit [51]) are valid in the region of the G(1.5-1.6) 
cannot be applied. 

4.2. Evidences that the a is a gluonium 

In section [21 we have given several indications that the 
a can be a gluonium from the values of its small 77 and 
large hadronic widths and of its on-shcll mass. In the 
previous section, the value of its K + K~ coupling, which 
is about 0.7 times of its tt + tt~ one (see Table [3]), gives 
a further support on its gluonium structure due to the 
fact that a tttt molecule and a tetraquark assignements 
would lead to a null value of the K + K~ coupling. This 
gluonium nature of the a can indicate that the effective 
KKtttt four-vertex used in a kaon loop model for suc- 
cessfully explaining e.g. 4> — > 7 7T7r can be induced by a 
s-channel exchange of a gluonium state. Indeed, using 
a gluonium picture, one can predict correctly this ob- 
served radiative width [21718] . For better comparing the 
results obtained in the complex plane with the QSSR 
results obtained in the real axis, we introduce like in 
[3] the on-shell meson masses and hadronic widths [52] . 
where the amplitude is purely imaginary at the phase 
90°: 



ReV({M° s ) 2 ) = M° s w 0.92 GeV . 



(23) 



In the same way as for the mass, one can define an 
"on-shell width" [see Eqs. (fT2"|) and (Tl3|) ] evaluated at 



— (Mos\2 . 



(m; 



M o T o 



Im V 

-Re V 



r 



0.7 GeV 



(24) 



which are comparable with the Breit-Wigner mass and 
width [45153154] : 



M BW » T B w » 1 GeV. 

These values lead to the on-shell coupling 

OS 

UaK+K- 



6 GeV 



tK 



0.8 



Q 



(25) 



(26) 



These fitted values and the one predicted in Eq. © 
strongly indicate a large gluonium component in the a 
wave function. 

5. Summary and conclusions 

• We have extracted the <jtt + tt~ and aK + K~ couplings 
from tttt — > tttt/KK scatterings using different models 



9 Notice that the analysis in Ref. [7] also indicates that erg (in the 
real axis) having a mass below 750 MeV cannot be wide (< 200 
MeV) due to the sensitivity of the coupling to M a . 
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of K- and S-matrices. We have shown in Table [3] the 
different results and have compared them with the ones 
from 4> — * 0"//o(98O) 7 and J/ip — ► <j> tttt/KK decays. 
A comparison of these results with the predictions from 
LET®QSSR favours a large gluonium/glueball content 
in the a wave function. 

• This feature can explain the splitting of its mass from 
the pion one, which is similar to the n—f]' mass splitting 
occuring in the well-known U(1)a anomaly 55J. The 
similarity with the rf is also signaled by the affinity of 
the a to couple to Goldstonc bosons indicating a large 
OZ1 violation in its decay into tttt and a large value of 
its KK coupling. In this way, the a associated to the 
trace of the energy- momentum tensor 9£ [U(l)v confor- 
mal anomaly] can be considered as the partner of the 77 
of the U(1)a anomaly. This non-gg property of the a 
may also go in lines with the \/N c counting done in [55] 
using unitarized ChPT partial waves for describing tttt 
scatterings. 

• Therefore, we consider that the assignement for the 
a as a non-strange partner of the k/ Kq(800) (if con- 
firmed), like often advocated in the four-quark litera- 
ture (see e.g. |57|lj ). may not be appropriate. The lat- 
ter having an isospin 1/2 cannot have a gluonium in its 
wave function. Indeed, the k can likely be the isoscalar 
partner of the ao(980), where, in the standard classifi- 
cation of current algebra, they are respectively the us 
and ud mesons having a mass around 1 GeV |8l7j which 
are associated to the divergence of the vector currents. 
Also, within this description, a successful determination 
of the running light quark mass differences from QSSR 
have been achieved [28|58|59j . 

• Finally, with the gluonium nature of the a and its 
large coupling to pseudoscalar boson pairs, the effective 
four-meson vertex KKitTt used e.g. for explaining the 
<\> — > 77T7T process can be due to a s-channel exchange of 
a glue rich state. We plan to come back to this point as 
well as to the study of the gluonium-quarkonium mixing 
below 1 GeV which has been initiated in |12l7j . 
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